Transport anomaly at the ordering transition for adatoms on graphene 
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We analyze a manifestation of the partial ordering transition of adatoms on graphene in resistiv- 
ity measurements. We find that Kekule mosaic ordering of adatoms increases sheet resistance of 
graphene, due to a gap opening in its spectrum, and that critical fluctuations of the order parameter 
lead to a non-monotonic temperature dependence of resistivity, with a cusp-like minimum at T = T c . 



Impurities in metals experience a long-range RKKY 
interaction due to polarization of the electron Fermi sea 
(Friedel oscillations) pQ. For surface adsorbents such an 
interaction may result in their structural ordering, re- 
peating the pattern of the Friedel oscillations of electron 
density [2J. In particular, a dilute ensemble of adatoms 
on graphene may undergo a partial ordering transition 

[3HB]. 

Unlike other materials, the RKKY interaction between 
adatoms on graphene exists even at zero carrier density, 
with a characteristic long-range 1/r 3 dependence, and it 
exhibits the Friedel oscillations which are commensurate 
with the underlying honeycomb lattice. For adatoms re- 
siding above the centers of the honeycombs, the interval- 
ley scattering of the electrons by adatoms leads to Friedel 
oscillations that resemble a y/3 x y/3 charge-density wave 
superlattice. Positions of each individual adatom rela- 
tive to such superlattice can be characterized by one of 
three vectors, u = (cos , sin ) , with m= —1,0, 1. 
The transition of an ensemble of adatoms into a "Kekule 
mosaic" ordered state [3] , characterized by the order pa- 
rameter u, falls in the symmetry class of 3-state Potts 
models [7]. 

In this paper we analyze how partial "Kekule" order- 
ing of adatoms on graphene affects its resistivity, p, in 
the regime of low coverage n^a 2 <C 1 (tii is the con- 
centration of adatoms, a is the lattice constant). The 
behaviour of the temperature dependent resistivity cor- 
rection 5p{T) = p(T) — p(oo) is sketched in Fig. [I] At 
T <C T c (region I) the temperature dependence of the re- 
sistivity is dominated by a non- vanishing order parameter 
u causing the amplified intervalley mixing and opening 
a gap A oc u ~ (T c — T) 13 in the corners of the Bril- 
louin zone. As the temperature increases from T = 
to T = T c the resistivity correction monotonically de- 
creases as (T c - T) 2f} . At T > T c , critical fluctuations 
of the order parameter, characterized by the correlation 
length £ cx \T — T c \~ v , preceding formation of the or- 
dered phase lead to a non-monotonic feature in 8p{T). 
At high temperature T > T c (region III), the construc- 
tive interference of electron waves, scattered by adatoms 
within ordered clusters of size £, enhances resistivity. The 
effect, which becomes stronger upon approaching T c , is 
similar to the critical opalescence [5] in materials under- 
going structural phase transition or resistivity anomaly 
in bulk metals with magnetic impurities undergoing fer- 
romagnetic transition [3]. This enhancement saturates 
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FIG. 1: The predicted anomaly in the temperature-dependent 
resistivity of graphene decorated with adatoms in the vicin- 
ity of the Kekule ordering transition. The inset illustrates 
the Kekule mosaic ordered state and the assignment of Potts 
"spin" m = —1,0, 1 to various hexagons in the \/3 x \/3 su- 
perlattice. 



when £ becomes comparable to the electron wavelength 
X F , X F rj £. In the region II of temperatures T — > T c + 0, 
where £ 3> Xf, scattering of electrons is affected only by 
the gradient of the fluctuating order parameter u. The 
resistivity is thus reduced and a cusp-shape minimum at 
T = T c should be expected. 

In the following we assume that the electron concen- 
tration n e — Air/Xp is not high, n e -C rij, but the elec- 
tron Fermi wavelength is shorter than its mean free path, 
Xf <C I. This assumption also implies that in the ordered 
phase ksT, A <C ep- 

The electrons are described by a four-component 
Dirac-like spinor * = [^Jk,A^K,B, ^K',b,^K',a] with 
the components corresponding to different valleys (K, 
K') and sublattices (A, B) [10]. In the absence of 
adatoms, quasiparticles are characterized by the linear 
spectrum |e p | = hvp and plane wave states (for e p > 0) 
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Here S is the area of the graphene sheet and p = 
(pcos(/?p,psin(/3 p ) is the electron wave vector. 
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The Hamiltonian describing graphene covered by 
adatoms has the form [3l [UJ 

H = ?w(pl!) + E7(r) + V"(r); (1) 
U = £ Itc(r - r,); V = RAwoE S,(uiA)J(r - r,). 

We assume that the dimensionless impurity potential A 
is small (A < 1) and will treat the electron-adatom in- 
teraction perturbatively. We use the set [5J |TT] of 4 x 4 
"sublattice" and "valley" matrices ^ x ,y,z and A-x,y,z 



A, 



) (T Z , Ay 



n, 



) tr,, A, 



n 



1 CT Z , 
> 00, 



where cr, and H, are Pauli matrices in sublattice (AB) 
and valley {KK') spaces. 

The form of the electron-adatom interaction in Eq. ([!]) 
is prescribed by the highly symmetric position of adatoms 
at the centers of hexagons. The [/-term does not violate 
the AB sublattice symmetry and scatters electrons with- 
out changing their valley state. The V^-term [2-4 is re- 
sponsible for the intervalley scattering of electrons. The 
sensitivity of the scattering phase of an electron to the 
position of the adatom in the y/3 x y/3 superlattice man- 
ifests itself by the Potts parameter u in the intervalley 
term. 

Each adatom creates Friedel oscillations of the elec- 
tron density leading to the RKKY-type interaction be- 
tween adatoms. The contribution of the [/-term, Eq. [l] 
to such an interaction is a 1 /r 3 repulsion independent on 
the Potts spins. At the same time the symmetry-breaking 
coupling leads to 3-state Potts model with a long-range 
interaction —Uj ■ u;/r? ; [3]. 

To minimize the interaction energy adatoms have to 
take equivalent positions within the V3 x y3 superlat- 
tice unit cells. A Monte Carlo simulation of the corre- 
sponding 1 /r 3 random-bond Potts model [3] has revealed 
such an ordered phase below the transition temperature 
T c « Q.Q\ 2 { ni a 2 f/ 2 hv/a. 

The effects of adatoms ordering on the electron trans- 
port are encoded in the correlation function, 



u a {v)uP{Y') - u a (r) u^(r') = S a/3 g(\r - r'|), 

for which the theory of critical phenomena predicts the 
scaling form [T^] 

g(r) = ^j=^ , K(y) = K 1 (y) + y^rK 2 (y). (2) 



Here, £ ~ n l 1/2 |(T -T c )/T c \~ u is the correlation length, 
and rj and v are critical exponents. It has been shown 
[7J that for the 3-state Potts model on a square lattice 
■q = 4/15, however, the value of r\ for random-bond Potts 
models with a long-range interaction is still unknown 
p~3|. IT4"] . In the critical region r -C £, the correlation 
function has essentially r~ v behaviour, with a correc- 
tion (second term) related to the specific heat anomaly 



C ~ \T-T c \~ a [7J. At large distances, r > £, g(r) decays 
exponentially, according to the Ornstein-Zernike theory 
[121 US] - Overall, for T > T c the scaling functions Ki^(y) 
have the following asymptotics: 

Kx(y < 1) « ci, K 2 (y<l)«-c 2 (ci,c 2 ~l); 

«i(y> !) « ^T/2=V «2(2/>l)~e v, (3) 

whereas for T < T c , 

«l(2/<l)~ci, K 2 (y < 1) ~ c 2 ; 
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At T < T c , the order parameter also acquires a ho- 
mogenious average u oc (T c — T) 13 . As long as u exceeds 
fluctuations of u (i.e. far enough from T c ), the elec- 
tron states, which wavelength is larger than the distance 

— 1/2 

between adatoms n i , are described by the effective 
mean-field Hamiltonian, 

H = hvpY^ + riihXva'S z (uA). 

Accordingly the spectrum e 2 , — (hvp) 2 + A 2 acquires a 
gap, 



A(T) w mh\va(l - T/T C )P, 



(4) 



such that A(0) ^> T c . The plane wave eigenstates of H 
are mixed between the two valleys and take the form (for 
£ P > 0) 



±1,P> = 



e D ±A 
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Intravalley and intervalley scattering determined by 
U(r) and F(r) in Eq. ^ respectively do not interfere 
with each other. Hence, the total momentum relaxation 
rate is the sum of the two electron scattering rates, 

T-^r^+rri, (5) 

where To and r, stand for intravalley and interval- 
ley momentum relaxation times. For the temperature- 
dependent Drude resistivity of graphene sheet (recall that 
ksT, A -C ep) we thus have 



(6) 



where vp — Hv 2 pf/£f is the Fermi velocity, v = 
2eF/{'Kh 2 v 2 ) is the density of states, and the Fermi en- 
ergy and momentum are related to the electron density 
as pf = \ZTTTie and Ef = \/irh 2 v 2 n e + A 2 . 
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The temperature dependence p(T) at T < T c is domi- 
nated by the effect of the order parameter u on the chiral 
plane wave functions and thus on the scattering rates, in 
particular Tq 1 . In the Born approximation 



approximate g(k) ~ g(0) oc (y^rii^) 2 11 and find that 
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where u>(fc) = J dv e w(r) and 99 is the scattering angle. 
The form- factor ro(tp) arises from the overlap integral 
between plane wave states and reflects the absence of the 
backscattering for A = 0. Thus, for T < T c , we find 



6p(T) _ 4A 2 (T) J** dtp w 2 (2p F sin sin 



2 '£ 
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p(oo) Tm e h 2 v 2 J 27r dip w 2 (2p F sin f ) sin 2 tp ' 



(8) 



The temperature dependence p(T) at T > T c is deter- 
mined by the effect of the ordering of adatoms on the 
intervalley scattering. Consider the scattering amplitude 



(K'p'\V\Kp) = 



hXva • fp 
2iS blU 



(9) 
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At temperatures far from T c , T 3> T c , adatom posi- 
tions on the superlattice are random so that mi take 
values —1,0 and 1 with equal probabilities. As a re- 
sult, the absolute value of scattering amplitude can be 
estimated as |( J K 7 p'|V'|.Kp)| ~ \friiX F . Upon approach- 
ing T c from above, clusters of ordered adatoms with a 

— 1/2 

characteristic size £ 3> n i start appearing. In the 
sum ([9]) such a cluster generates constructive interfer- 
ence between terms with the same value of mi provided 
that £ < Xf- This increases the scattering amplitude, 
I (-fC'p'lV"! A'p) I ~ n^Xp. A further increase of the corre- 
lation length, £ > Xf, has an opposite effect on scatter- 
ing: electrons get scattered only by the gradients in the 
smoothly fluctuating field u. 

The intervalley momentum relaxation rate (both at 
T < T c and T > T c ) [T7] t" 1 can be expressed in 
terms of the Fourier transform of the correlation func- 
tion, g{k) = J dre tkr g(r), 
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At T > T c (A = 0) temperature dependence p(T) 
comes from the correlation function g(k) in Eq. ( 10 ). Far 



from the phase transition, \T — T c \ ~ T c , where £ < Xp 
(region III in Fig. [I]), electrons are effectively scattered 
by small clusters of ordered adatoms. In this region we 



where C = (3tt 2 /2) £ dy y 1 ' 11 K (y) is a dimensionless 
constant. 

In the vicinity of the critical point, such that Xf < £ 
(region II in Fig. [I]), electrons experience multiple scat- 
terings within one cluster with a small wave vector trans- 
fer, ~ This makes p(T) sensitive to the critical be- 
haviour of the correlation function at r < £. This region 
is easie r to analyze by performing the angular integration 
in Eq. (10) and expressing t" 1 in terms of the function 
n(y) denned in Eq. 
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Here, / can be expressed in terms of Bessel functions as 
f(x) = -[x 2 J$(x) - x 2 Jl{x) + J 2 (x) - xJ a (x)J 1 (x)]. 



To evaluate the integral in Eq. ( 12 ) we divide the inte- 



gration interval [0, 00] into two parts, [0, j/o] an( i [z/o? °°] ; 
where 1 » yo ^> 1/pf£, — > 0. For the interval [j/0,00], we 
use the fact that fipF^y) is a fast oscillating function, 
f(x > 1) w 2sin(2.T) + 3 cos(2a;)/(2x), and that 
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F(y) sin(Ay)dy = 



F(y ) cos(Ay Q ) 
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for A ^> 1 and F(oo) = 0. For the interval [0, yo] m 
the leading order in 1/Vf£, the result is determined by 
the values of /«i(0) and «2(0) in Eq. For this, we 
expand K\ and K2 (which vary at the scale of y ~ 1) 
into Taylor series, evaluate the corresponding integrals 
in the leading orders in yo ^ 1) & n d combine with the 
contribution from the interval [2/0,00]- As a result, the 
term with K\ in Eq. (12 1 produces a finite contribution 



when pf£, —> 00, and we find that 

c^lA 2 (0) 

p(T c ) - p(oo) = B(ri) 

e z v z hni 
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(13) 



where B(x) = xT (3±») T (l - §) /T (l + f ) T (2 + f). 
The next term in the expansion k\ (y) — K\ (0) + yn'i (0) + 
... generates a contribution 0{l/pp^), which, for T — > 
T c , is less relevant than a more singular (T — In- 
dependent |14| contribution from the K2 term in Eq. ( |12[ ) . 
Following to same steps, we find that the latter term gives 
rise to the cusp in the p{T) dependence in the region II 
near T c in Fig. [l] 



p(T)-p(T c ) = c 2 B( V -2 1 ) 
p(T c )-p(oo) Cl B(r/)(™ e £ 2 )T 



cx (T-T c 



\l-ct 



(14) 
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where 7 = (1 — a)/2v. 

The behaviour of 5p(T) at T < T c (region I in Fig. [if 
is determined by two contributions. One part, Spi/p oc 
(T c — T) 1_Q , is related to the specific heat anomaly cor- 
rection to the correlation function and can be obtained 
the same way as Eq. (14 1. The other contribution, 
T) 2/3 , is due to the formation of a non- 



5p 2 /p oc (T c 

zero order parameter in the Kekule-ordered phase. The 
second correction dominates when 2/3 < 1 — a, which is 
the case for the expected values of critical exponents [H] . 
As a result, we attribute the rise of resistivity at T < T c 
near the cusp at T = T c to the formation of a spectral 
gap in graphene due to the Kekule mosaic ordering. The 
resulting behaviour of resistivity as a function of temper- 
ature for all three regimes is plotted in Fig. [I for C\ = 0.5, 
c 2 = 0.15, e F = QAv^/W u ni \ 2 a 2 = 0.005 [H]. 

In conclusion, we investigated electron transport in 
graphene covered by a dilute ensemble of adatoms re- 



siding over the centers of hexagons. We calculated the 
temperature dependence of the resistivity p{T), which 
appears to be non-monotonic and has a non-analytic cusp 
at T = T c . Since the form of the cusp depends on the crit- 
ical indices a and f3 of the phase transition, experimental 
observation of such anomaly may facilitate their measure- 
ments. The form of p(T) shown in Fig. [I] appears to be 
generic for partially ordered dilute ensembles of adatoms 
with alternative positioning on the honeycomb lattice: 
(i) over the sites O [6] and (ii) over carbon-carbon bonds 
[I]. Since those also fall into the class of Potts models 
[(i) - 2 value Potts model and (ii) - 3 value Potts model], 
the anomalies in p(T) can be described using Eqs. fl8|i"4[ ) 
with appropriate critical indices, a and /3. 
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